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Abstract. We show the existence of a fully polynomial-time approximation 
scheme (FPTAS) for the problem of maximizing a non-negative polynomial 
over mixed- integer sets in convex polytopcs, when the number of variables 
is fixed. Moreover, using a weaker notion of approximation, we show the 
existence of a fully polynomial-time approximation scheme for the problem of 
maximizing or minimizing an arbitrary polynomial over mixed-integer sets in 
convex polytopes, when the number of variables is fixed. 



1. Introduction 

A well-known result by H.W. Lenstra Jr. states that integer linear programming 
problems with a fixed number of variables can be solved in polynomial time on 
the input size |12j . Likewise, mixed integer linear programming problems with a 
fixed number of integer variables can be solved in polynomial time. It is a natural 
question to ask what is the computational complexity, when the number of variables 
(or the number of integer variables) is fixed, of the non-lin ear mixed integer problem 

max f(xt, ... ,x dl ,Zi, ... ,z d2 ) (la) 

s.t. Ax + Bz<h (lb) 

ij£R for i = 1, . . . , d\, (lc) 

Zi G Z for i = 1, . . . , d 2 , (Id) 

where / is a polynomial function of maximum total degree D with rational coeffi- 
cients, and A £ Z pxdl , B G Z pxd2 , b G Z p . We are interested in general polynomial 
objective functions / without any convexity assumptions. 
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Throughout the paper we assume that the inequality system An + Bz < b 
describes a convex polytope, i.e., a bounded polyhedron, which we denote by P. The 
reason for this restriction are fundamental noncomputability results for problems 
involving polynomials and integer variables. Indeed, when we permit unbounded 
feasible regions, there cannot exist any algorithm to decide whether there exists 
a feasible solution to |T]) with /(x, z) > a (for a prescribed bound a), ruling out 
the existence of an optimization algorithm or any approximation scheme. This is 
due to the negative answer to Hilbert's tenth problem by Matiyasevich [T51 114|. 
Due to Jones' strengthening of this negative result [11) . there also cannot exist any 
such algorithm for the cases of unbounded feasible regions for any fixed number of 
integer variables g?2 > 10. 

For the purpose of complexity analysis, we assume that the data A, B, and b are 
given by the binary encoding scheme, and that the objective function / is given as 
a list of monomials, where the coefficients are encoded using the binary encoding 
scheme and the exponent vectors are encoded using the unary encoding scheme. In 
other words, the running times arc permitted to grow polynomially not only in the 
binary encoding of all the problem data, but also in the maximum total degree D 
of the objective function /. 

It is well-known that pure continuous polynomial optimization over polytopes 
(d-2 = 0) in varying dimension is NP-hard and that a fully polynomial time ap- 
proximation scheme (FPTAS) is not possible (unless P = NP). Indeed the max-cut 
problem can be modeled as minimizing a quadratic form over the cube [— 1,1] . 
Hastad [10) proved that the max-cut problem cannot be approximated to a ratio 
better than 1.0625 (unless P = NP). This excludes the possibility of a polynomial 
time approximation scheme for JTJ in varying dimension, even when the number of 
integer variables is fixed. 

On the other hand, pure continuous polynomial optimization problems over poly- 
topes (o?2 = 0) can be solved in polynomial time when the dimension d\ is fixed. 
This follows from a much more general result on the computational complexity of 
approximating the solutions to general algebraic formulae over the real numbers by 
Renegar [19]; see also [TBI [T7| [T8] . 

However, when we permit integer variables (c?2 > 0), it turns out that, even for 
fixed dimension d\ + d-2 = 2 and objective functions / of maximum total degree 
D = 4, problem {TJ is an NP-hard problem [7j. Thus the best we can hope for, 
even when the number of both the continuous and the integer variables is fixed, is 
an approximation result. This paper presents the best possible such result: 

Theorem 1 (Fully polynomial-time approximation scheme). Let the dimension 
d = d\ + g?2 be fixed. 

(a) There exists a fully polynomial time approximation scheme (FPTAS) for the 
optimization problem ([JJ for all polynomial functions f(x\, . . . , Xd x ,Z\,..., Zd 2 ) 
with rational coefficients that are non-negative on the feasible region (jlbHld[) . 

(b) Moreover, the restriction to non-negative polynomials is necessary, as there does 
not even exist a polynomial time approximation scheme (PTAS) for the maxi- 
mization of arbitrary polynomials over mixed-integer sets in polytopes, even for 
fixed dimension d > 2, unless P = NP. 

The proof of IThcorcm II is presented in lscction 51 As wc will sec. IThcorcm H is 
a non-trivial consequence of the existence of FPTAS for the problem of maximizing 
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a non-negative polynomial with integer coefficients over the lattice points of a con- 
vex rational polytope. That such FPTAS indeed exist was recently settled in our 
paper [7J. The knowledge of paper [7J is not necessary to understand this paper 
but, for convenience of the reader, we include a short summary in Isection 21 Our 
arguments, however, are independent of which FPTAS is used in the integral case. 

Our main approach is to use grid refinement in order to approximate the mixed- 
integer optimal value via auxiliary pure integer problems. One of the difficulties 
on constructing approximations is the fact that not every sequence of grids whose 
widths converge to zero leads to a convergent sequence of optimal solutions of grid 
optimization problems. This difficulty is addressed in Isection 31 In Isection 4l we 
develop techniques for bounding differences of polynomial function values. Section[5] 
contains the proof of lThcorcm 11 

Finally, in Isection 6l we study a different notion of approximation. The usual 
definition of an FPTAS uses the notion of e-approximation that is common when 
considering combinatorial optimization problems, where the approximation error is 
compared to the optimal solution value, 

|/(^e; 2 e ) f (x max , Z max ) | € e /'(x max , Z max ), (2) 

where (x e ,z e ) denotes an approximate solution and (x max ,z max ) denotes a maxi- 
mizer of the objective function. In Isection 6l we now compare the approximation 
error to the range of the objective function on the feasible region, 

| f(x e ,z e ) f (x max , z max ) I 6 1 /*(x max , Z max ) ,/(x m i n , Z m j n ) I , (3) 

where additionally (x m i n ,z m i n ) denotes a minimizer of the objective function on 
the feasible region. This notion of approximation was proposed by various authors 
[20l [3l [5j. It enables us to study objective functions that are not restricted to be 
non-negative on the feasible region. We remark that, when the objective function 
can take negative values on the feasible region, ([3]) is weaker than ([2]). Therefore 
IThcorcm ll fb) does not rule out the existence of an FPTAS with respect to this 
notion of approximation. Indeed we prove: 

Theorem 2 (Fully polynomial-time weak- approximation scheme). Let the dimen- 
sion d = d\ + di be fixed. Let f be an arbitrary polynomial function with rational 
coefficients and maximum total degree D, and let P C R d be a rational convex 
polytope. 

(a) In time polynomial in the input size and D, it is possible to decide whether f 
is constant on P n (R dl xZ*). 

(b) In time polynomial in the input size, D, and - it is possible to compute a 
solution (x £ ,z e ) £ Ffl (R dl x Z* 2 ) with 

| f (-X-e ; ^e) / (^max s ^max) | — ^ | f (^maxi ^max) f (-^min; ^min) | ■ 

Notation. As usual, we denote by Q[xi, . . . , Xd t , z±, . . . , Zd 2 ] the ring of multivariate 
polynomials with rational coefficients. For writing multivariate polynomials, we 
frequently use the multi- exponent notation, z a = z" 1 ■ ■ ■ z^ d . 

2. An FPTAS for the integer case 

The first fully polynomial-time approximation scheme for the integer case ap- 
peared in our paper [7J. It is based on Alexander Barvinok's theory for encoding 
all the lattice points of a polyhedron in terms of short rational functions [TJ [5] . The 
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set P n Z d is represented by a Laurent polynomial gp{z) = X)aePnz d zQ - From 
Barvinok's theory this exponentially-large sum of monomials gp(z) can be written 
as a polynomial-size sum of rational functions (assuming the dimension d is fixed) 
of the form: 

where / is a polynomial-size indexing set, and where Ei £ {1,-1} and U;, Vjj £ Z d 
for all i and j . There is a polynomial-time algorithm for computing this represen- 
tation m |3 El El . 

By symbolically applying differential operators to the representation ((4]) . we can 
compute a short rational function representation of the Laurent polynomial 

g PJ (z)= J2 /(") z " ( 5 ) 
aePnz d 

In fixed dimension, the size of the expressions occuring in the symbolic calculation 
can be bounded polynomially: 

Lemma 3 ([7], Lemma 3.1). Let the dimension d be fixed. Let gp(z) = ^2 a£PnZd z a 
be the Barvinok representation of the generating function of P Z d . Let f £ 
Zi[xi, . . . ,Xd] be a polynomial of maximum total degree D. We can compute, in 
time polynomial in D and the input size, a Barvinok representation gpj(z) for the 
generating function X) a ePnz d f( a ) za - 

Now we present the algorithm to obtain bounds Uk, Lfe that reach the optimum. 
We make use of the elementary fact that, for a set S = {si, . . . , s r } of non-negative 
real numbers, 

max{si, . . . , s r } = lim \J s\ + ■ ■ ■ + s r k . (6) 

k — >oo 

Algorithm 4 (Computation of bounds). 

Input: A rational convex polytopc P C R d , a polynomial objective / £ Z[xi, . . . , Xd) 
of maximum total degree D that is non- negative over P n Z d . 
Output: A nondecreasing sequence of lower bounds L^, and a nonincreasing se- 
quence of upper bounds Uk , both reaching the maximal function value /* of / over 
P fl Z d in a finite number of steps. 

(1) Compute a short rational function expression for the generating function 
gp(z) = EaePnz^ z "- Using residue techniques, compute |PnZ d | = gp(l) 
from gp(z). 

(2) From the rational function gp(z) compute the rational function represen- 
tation of g P jk(z) of J2aePnz d f k ( a ) za bv lLcmma 31 Using residue tech- 
niques, compute 



L 



k • = 



y/9p,f'(l)/9pj°(l) and U k ~ ^5^(1) J 



Theorem 5 ([7], Lemma 3.3 and Theorem 1.1). Let the dimension d be fixed. Let 
P C R d be a rational convex polytope. Let f be a polynomial with integer coefficients 
and maximum total degree D that is non-negative on P P\Z d . 

(i) \Algorithm 4\ computes the bounds Lk, Uk in time polynomial in k, the input 
size of P and f, and the total degree D. The bounds satisfy the following 
inequality: 



u k -L k <r- {\Pnz d \-i 
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(ii) For k = (1 + 1/e) log(|P (~l Z d |) ('a number bounded by a polynomial in the 
input size), Lk is a (1 — ^-approximation to the optimal value f* and it can 
be computed in time polynomial in the input size, the total degree D, and 1/e. 
Similarly, Uk gives a (1 + ^-approximation to /*. 

(iii) With the same complexity, by iterated bisection of P, we can also find a feasible 
solution x £ G P H Z d with 

|/(x e )-/*| <ef*. 

3. Grid approximation results 

An important step in the development of an FPTAS for the mixed-integer opti- 
mization problem is the reduction of the mixed-integer problem ([T]) to an auxiliary 
optimization problem over a lattice -^Z dl x Z d2 . To this end, we consider the grid 
problem with grid size to, 

max f(x 1 ,...,x dl ,Zi,...,Zci 2 ) 

s.t. Ax + Bz < b 

i (7) 
x t G — Z for i = 1,.. .,d\, 

Zi G Z for % = 1 , . . . , d% ■ 

We can solve this problem approximately using the integer FPTAS (jThcorcm 51) : 

Corollary 6. For fixed dimension d — d± + di there exists an algorithm with 
running time polynomial in log to, the encoding length of f and of P, the maximum 
total degree D of f, and ^ for computing a feasible solution (x™,z™) G Pn(^-Z dl x 
Z rf2 ) to the grid problem ([7]) with an objective function f that is non-negative on 
the feasible region, with 

j-(x™zr)>(l-c)/(x"\z™), (8) 
where (x m ,z m ) G Pf) (^Z dl xZ 4 ) is an optimal solution to 0. 
Proof. We apply [Theorem 5l to the pure integer optimization problem: 
max /(x, z) 

s.t. Ax + mBz, < mh , , 

(9) 

ii G Z for i = 1 , . . . , di , 

Zi G Z for i = 1 , . . . , g?2 7 

where /(x,z) := m D f(^Si,z) is a polynomial function with integer coefficients. 
Clearly the binary encoding length of the coefficients of / increases by at most 
[I? log to], compared to the coefficients of /. Likewise, the encoding length of 
the coefficients of mB and mh increases by at most [log to]. By Theorem 1.1 
of [7], there exists an algorithm with running time polynomial in the encoding 
length of / and of Ax + m_Bz < mh, the maximum total degree D, and - for 
computing a feasible solution (x™, z™) G Pfl (^Z dl x Z da ) such that /(x™, zf ) > 
(1 - e)/(x m , z m ), which implies the estimate ©. □ 

One might be tempted to think that for large-enough choice of to, we immediately 
obtain an approximation to the mixed-integer optimum with arbitrary precision. 
However, this is not true, as the following example demonstrates. 
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Figure 1 . A sequence of optimal solutions to grid problems with 
two limit points, for even m and for odd m 

Example 7. Consider the mixed-integer optimization problem 

max 2z — x 
s.t. z < 2x 

z<2(l-x) < 10 > 
x G R> , z G {0, 1}, 

whose feasible region consists of the point (^,1) and the segment {(x,0) : x G 
[0,1]}. The unique optimal solution to ([10]) 1. Now consider the 

sequence of grid approximations of (fT0|l where x G ^Z>o- For even m, the unique 
optimal solution to the grid approximation is x = |, z = 1. However, for odd m, 
the unique optimal solution is x = 0, z = 0. Thus the full sequence of the optimal 
solutions to the grid approximations does not converge since it has two limit points; 
see |Figure 1| 

Even though taking the limit does not work, taking the upper limit does. More 
strongly, we can prove that it is possible to construct, in polynomial time, a sub- 
sequence of finer and finer grids that contain a lattice point (|_x*](5,z*) that is 
arbitrarily close to the mixed-integer optimum (x*,z*). This is the central state- 
ment of this section and a basic building block of the approximation result. 

Theorem 8 (Grid Approximation). Let d\ be fixed. Let P = { (x, z) G R dl+d2 : 
Ax + Bz < b}, where A G Z' pxdl , B G Z pxd2 . Let M G R be given such that 
P Q { ( x , z ) G R dl + d2 : \xi\ < M for i — 1, ...,d\}. There exists a polynomial- 
time algorithm to compute a number A such that for every (x*, z*) G Pl~l(R dl x Z d2 ) 
and 5 > the following property holds: 

Every lattice ^;Z dl for m = kA and k > %d\M contains a lattice 
point [x.*1s such that ([x*]s,z*) G Pn(±Z dl xZ d2 ) and ||lx*]<5- 
x* II < S. 

II oo — 

The geometry of IThcorcm 81 is illustrated in |Figurc 2| The notation [ x *]<5 nas 
been chosen to suggest that the coordinates of x* have been "rounded" to obtain 
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Figure 2. The principle of grid approximation. Since we can 
refine the grid only in the direction of the continuous variables, we 
need to construct an approximating grid point (x, z* ) in the same 
integral slice as the target point (x*, z*). 

a nearby lattice point. The rounding method is provided by the next two lemmas; 
IThcorcm 81 follows directly from them. 

Lemma 9 (Integral Scaling Lemma). Let P = { (x, z) G R d i+*i ; Ax + Bz <b}, 
where A G Z pxdl , B G Z p 2 . For fixed di, there exists a polynomial time algorithm 
to compute a number A G Z>q such that for every z G Z d2 the polytope 



is integral, i.e., all vertices have integer coordinates. In particular, the number A 
has an encoding length that is bounded by a polynomial in the encoding length of P. 

Proof. Because the dimension d\ is fixed, there exist only polynomially many sim- 
plex bases of the inequality system Ax < b — Bz, and they can be enumerated 
in polynomial time. The determinant of each simplex basis can be computed in 
polynomial time. Then A can be chosen as the least common multiple of all these 
determinants. □ 

Lemma 10. Let Q C R d be an integral polytope. Let M G R. be such that Q C 
{ x G R d : \xi\ < M for i = 1, . .. ,d}. Let x* G Q and let S > 0. Then every lattice 
j:7i d for k > jdM contains a lattice point x G Q H \'L d with ||x — x*||oo < 5. 

Proof. By Carathcodory's Theorem, there exist d+1 vertices x°, . . . , x d G Z d of Q 
and convex multipliers Aq, . . . , such that x* = X^=o ^i xI - Let A^ := v [k\\ > 
for i = 1, . . . , d and Ag := 1 — $^ =1 A£ > 0. Moreover, we conclude A; — X[ < i for 



i = 1, . . . , d and A - A = £ti( A * - K) < 4- Thcn x : = Etc e n 



AP z = {Ax: (x,z) 6P} 



and we have 



x 



x 



< 



i=0 
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which proves the lemma. 



□ 



4. Bounding techniques 

Using the results of lsection 31 we are now able to approximate the mixed-integer 
optimal point by a point of a suitably fine lattice. The question arises how we can 
use the geometric distance of these two points to estimate the difference in objec- 
tive function values. We prove [Lemma 1 H that provides us with a local Lipschitz 
constant for the polynomial to be maximized. 

Lemma 11 (Local Lipschitz constant). Let f be a polynomial in d variables with 
maximum total degree D. Let C denote the largest absolute value of a coefficient 
of f . Then there exists a Lipschitz constant L such that |/(x) — /(y)| < L||x — y||oo 
for all \xi\, \yi\ < M. The constant L is 0{D d+l CM D ). 

Proof. Let /(x) = J2 a ev CaX - a > where T> C Z> is the set of exponent vectors of 
monomials appearing in /. Let r = \T>\ be the number of monomials of /. Then 
we have 

|/(x)-/(y)|<£|c Q ||x a -y a |. 

ct#0 

We estimate all summands separately. Let a ^ be an exponent vector with 
n := Y^i=i a i ^ D- Let 

a = cx° > a 1 > ■ ■ ■ > ot n = 



be a decreasing chain of exponent vectors with a 1 



i Ji for i 



1. 



Let [3 l 
sum" 



a — a 1 for i 



E 

i=l 
n 

E 

i=l 



= o, 



, n. Then x a — y a can be expressed as the "telescope 



rP^ 



rP' 



Since |x a 'y' 3 * 



< M n 1 and n < D, we obtain 
|x a -y a | <£>.||x-y|| 



thus 



rD-l\ 



Let L := CrDM 



D-l 



|/(x)-/(y)| KCrDM 1 - -||x-y| 
Now, since r = 0(D d ), we have L = 0(D d+1 CM D ). 



□ 



Moreover, in order to obtain an FPTAS, we need to put differences of function 
values in relation to the maximum function value. To do this, we need to deal 
with the special case of polynomials that are constant on the feasible region; here 
trivially every feasible solution is optimal. For non-constant polynomials, we can 
prove a lower bound on the maximum function value. The technique is to bound the 
difference of the minimum and the maximum function value on the mixed-integer 
set from below; if the polynomial is non-constant, this implies, for a non-negative 
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Figure 3. The geometry of ILcmma 131 For a polynomial with 
maximum total degree of 2, we construct a refinement j,Z d (small 
circles) of the standard lattice (large circles) such that P fl \Z d 
contains an affine image of the set {0, 1, 2} d (large dots). 



polynomial, a lower bound on the maximum function value. We will need a simple 
fact about the roots of multivariate polynomials. 

Lemma 12. Let f G Q[xi, . . . ,Xd] be a polynomial and let D be the largest power 
of any variable that appears in f. Then f = if and only if f vanishes on the set 
{0,...,D} d . 

Proof. This is a simple consequence of the Fundamental Theorem of Algebra. See, 
for instance, [H Chapter 1, §1, Exercise 6 b]. □ 

Lemma 13. Let f G Q[ii, ■ • • , xj] be a polynomial with maximum total degree D. 
Let Q C R d be an integral polytope of dimension d' < d. Let k > D d 1 . Then f is 
constant on Q if and only if f is constant on Q D \Z d . 

Proof. Let x° G Q(~)Z d be an arbitrary vertex of Q. There exist vertices x 1 , . . . , x d G 
Q n Z d such that the vectors x 1 — x°, . . . , x rf — x° G Z d arc linearly independent. 
By convexity, Q contains the parallelepiped 

S := { x° + Eti M* 1 - x°) : A,- G [0, for i = 1, . . . , d' }. 

We consider the set 

S k = \Z d n S D { x° + Eti f(x l - x°) : m G {0, 1, . . . ,D} for i = 1, . . . , df }; 

see |Figure 3| Now if there exists a c G R with /(x) = c for all x G Q fl \Z d , then 
all the points in Sk are roots of the polynomial / — c, which has only maximum 
total degree D. By ILcmma 121 (after an affine transformation), / — c is zero on the 
affine hull of Sk] hence / is constant on the polytope Q. □ 

Theorem 14. Let f G Z[x\, . . . , Xd xl z\, . . . , Zd 2 ]- Let P be a rational convex poly- 
tope, and let A be the number from [Lemma (A Let m = kA with k > Dd\, k G Z. 
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Optimal mixed-integer solution 
Optimal grid solution 
Rounded mixed-integer solution 

Approximative grid solution 



Figure 4. Estimates in the proof of lTheorem 11 (a) 

Then f is constant on the feasible region Pfl (R dl x Z d2 ) if and only if f is constant 
on Pn(±Z dl xZ d2 ). If f is not constant, then 

|/(x ma xi Z max ) /(x m i n , Z m i n ) | > HI , (H) 

where (x max , z max ) is an optimal solution to the maximization problem over the 
feasible region Pfl (R dl x Z d2 ) and (x m i n ,z m j n ) is an optimal solution to the min- 
imization problem. 

Proof. Let / be constant on PC) (^Z dl x Z d2 ). For fixed integer part z € Z d2 , we 
consider the polytope AP Z = { Ax : (x, z) G P}, which is a slice of P scaled to 
become an integral polytope. By applying ILcmma 131 with k = (D + l)d on every 
polytope AP Z , we obtain that / is constant on every slice P z . Because / is also 
constant on the set P (~l (^;Z dl x Z d2 ), which contains a point of every non-empty 
slice P z , it follows that / is constant on P. 

If / is not constant, there exist (x^z 1 ), (x 2 ,z 2 ) E P n (^Z dl x Z d2 ) with 
/(x^z 1 ) 7^ /(x 2 ,z 2 ). By the integrality of all coefficients of /, we obtain the 
estimate 

(/(xSz 1 )-/^ 2 )] >m- D . 

Because (x^z 1 ), (x 2 ,z 2 ) are both feasible solutions to the maximization problem 
and the minimization problem, this implies (jlip . □ 

5. Proof of ITheorem II 

Now we are in the position to prove the main result. 

Proof of \Theorem 1\ Part (a). Let (x*,z*) denote an optimal solution to the 
mixed-integer problem (JTJ) . Let e > 0. Wc show that, in time polynomial in 
the input length, the maximum total degree, and -, we can compute a point (x, z) 
that satisfies (jlbllldp such that 

|/(x,z)-/(x*,z*)|<e/(x*,z*). (12) 

We prove this by establishing several estimates, which are illustrated in |Figure~4) 
First we note that we can restrict ourselves to the case of polynomials with 
integer coefficients, simply by multiplying / with the least common multiple of all 
denominators of the coefficients. We next establish a lower bound on /(x*, z*). To 
this end, let A be the integer from ILcmma 9| which can be computed in polynomial 



/(x*,z*) 
/(x"\z m ) 



I 1 



<L5< f/(x*,z*) 
<§/(x"\z m )<f/(x*,z*) 



/(x e /2,Z E/2 
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time. Bv lTheorem 141 with m = DdiA, cither / is constant on the feasible region, 



or 



/(x*,**)>(DdiA)- D , (13) 
where D is the maximum total degree of /. Now let 

e 



' 2{Dd 1 A) D L(C,D,M) 
and let us choose the grid size 

~4 



A 



-(Z?di A) D L(C, D, M)d\M 

e 



(14) 



(15) 



where L{C,D,M) is the Lipschitz constant from [Lemma 1 ll Then we have m > 
AfdjM, so bv IThcorcm~8l there is a point ([x*] 5 ,z*) e Ffl (^Z dl x Z d2 ) with 
|| |_x*] g — x*|| < (5. Let (x m ,z m ) denote an optimal solution to the grid prob- 
lem (J7J). Because (|x*],5,z*) is a feasible solution to the grid problem (J7J), we have 

/(Lx*1j,z*) < /(x m ,z m ) </(x*,z*). (16) 

Now we can estimate 

|/(x*,z*)-/(x"\z™)| < |/(x*,z*)-/(Lx*l 5 ,z*)| 
<L{C,D,M)\\x* - Lx*l«L 
< L(C,D,M)8 

= e -{Dd 1 A)- D 

<f/(x*,z*), (17) 

where the last estimate is given by (|13fl in the case that / is not constant on the 
feasible region. On the other hand, if / is constant, the estimate (|17|l holds trivially. 

By Corollary 6 we can compute a point (x™ 2 ,z™ 2 ) ePn x Z d2 ) such 

that 

(1 - f)/(x™,z m ) < /(x™ /2 ,z^ 2 ) < /(x m ,z™) (18) 

in time polynomial in logm, the encoding length of / and P, the maximum total 
degree D, and 1/e. Here logm is bounded by a polynomial in log M, D and logC, 
so we can compute ( x ™2' z "/2) m time polynomial in the input size, the maximum 
total degree D, and 1/e. Now, using (fT5|) and (flT)) . we can estimate 

/( X *,z*) - /(xr /2 ,z:; 2 ) < /(x*, z *) - (i - §)/(x-\z m ) 

= f/(x*,z*) + (1 - !)(/(x*,z*) - /(x™,z™)) 
< f/(x*,z*) + f/(x%z*) 
= e/(x*,z*). 

Hence /(x^ 2 , z^ 2 ) > (1 - e)/(x*, z»). 

Part (b). Let the dimension d > 2 be fixed. We prove that there does not exist 
a PTAS for the maximization of arbitrary polynomials over mixed-integer sets of 
polytopes. We use the NP-complete problem AN1 on page 249 of [9]. This is to 
decide whether, given three positive integers a, 6, c, there exists a positive integer 
x < c such that x 2 = a (mod b). This problem is equivalent to asking whether the 



12 



J. A. DE LOERA, R. HEMMECKE, M. KOPPE, AND R. WEISMANTEL 



maximum of the quartic polynomial function f(x,y) = — (x 2 — a — by) 2 over the 
lattice points of the rectangle 

{1 — a (c — I) 2 — a 1 
(x,y) : 1 < x < c - 1, — — < y < ± ± \ 

is zero or not. If there existed a PTAS for the maximization of arbitrary polynomials 
over mixed- integer sets of polytopes, we could, for any fixed < e < 1, compute 
in polynomial time a solution (x t ,y e ) £ PflZ 2 with \f(x e ,y e ) — f(x*,y*)\ < 
e\f(x*, y*)\, where (x* , y*) denotes an optimal solution. Thus, we have f(x e , y t ) = 
if and only if f(x*,y*) = 0; this means we could solve the problem AN1 in 
polynomial time. □ 



6. Extension to arbitrary polynomials 

In this section wc drop the requirement of the polynomial being positive over the 
feasible region. As we showed in lThcorcm 11 there does not exist a PTAS for the 
maximization of an arbitrary polynomial over polytopes in fixed dimension. Wc 
will instead show an approximation result like the one in [5], i.e., we compute a 
solution (x e ,z € ) such that 

| f(x e ,Z e ) f (x max , Z max ) I € | ,/*(x max , Z max ) /(x m i n , Z m i n ) | , (19) 

where (x max ,z max ) is an optimal solution to the maximization problem over the 
feasible region and (x m ; n , z m ; n ) is an optimal solution to the minimization problem. 
Our algorithm has a running time that is polynomial in the input size, the maximum 
total degree of /, and -. This means that while the result of [5] was a weak version 
of a PTAS (for fixed degree), our result is a weak version of an FPTAS (for fixed 
dimension) . 

The approximation algorithms for the integer case (jThcorcm 5| and the mixed- 
integer case ([Theorem lj) only work for polynomial objective functions that are non- 
negative on the feasible region. In order to apply them to an arbitrary polynomial 
objective function /, we need to add a constant term to / that is large enough. As 
proposed in [7], we can use linear programming techniques to obtain a bound M 
on the variables and then estimate 

/(x) > -rCM D =: L , 

where C is the largest absolute value of a coefficient, r is the number of monomials 
of /, and D is the maximum total degree. However, the range |/(x max , z max ) — 
) can be exponentially small compared to L , so in order to obtain 
an approximation (x e ,z e ) satisfying (1 . we would need an (1 — e')-approximation 
to the problem of maximizing g(x, z) :~ /(x, z) — Lq with an exponentially small 
value of e'. 

To address this difficulty, we will first apply an algorithm which will compute 
an approximation [L h Ui\ of the range [/(x min , z min ), /( x max , z max )] with constant 
quality. To this end, we first prove a simple corollary of lTheorem 11 

Corollary 15 (Computation of upper bounds for mixed-integer problems). Let the 

dimension d — d\ + d,2 be fixed. Let P C R d be a rational convex polytope. Let 
f € Z[xi, . . . , Xd iy Z\, ... , Zd 2 ] be a polynomial function with integer coefficients and 
maximum total degree D that is non-negative on PC\ (R dl x Z d2 ) . Let 6 > 0. There 
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exists an algorithm with running time polynomial in the input size, D, and y for 
computing an upper bound u such that 

) <u< (l + 5)/(x max ,z max ), (20) 

where (x max , z max ) is an optimal solution to the maximization problem of f over 
PD (R dl x Z d2 ). 

Proof. Let e = By IThcorcm 11 we can, in time polynomial in the input size, 
D, and - = 1 + i, compute a solution (x e , z e ) with 

| /(^maxi Z max ) /(x € ,Z e )| < ^/(x max , Z max ). (^1) 

Let u := ^/(x, z e ) = (1 + 5)/(x e ,z e ). Then 



and 



/(x max ,z max ) < :-— /(x e ,z e ) = u (22) 



(! + *)/( x maxj z max) (1 + <^J /( x e 7 z e) 

= (l + <5)(l-e)u 

= (1 + 8)1 1- ]u = u. (23) 



1 

This proves the estimate (|20p . □ 
Algorithm 16 (Range approximation). 

Input: Mixed-integer polynomial optimization problem |T]), a number < 6 < 1. 
Output: Sequences {Li}, {[/;} of lower and upper bounds of / over the feasible 
region P n (R dl x Z d2 ) such that 

Li /(x m i n , Z m j n ) /(x max , z max ) < C/j (24) 

and 

lim - ^ | = c(/(x max ,z max ) - /(x min ,z min )), (25) 

i— >oo 

where c depends only on the choice of 8. 

(1) By solving 2d linear programs over P, we find lower and upper integer 
bounds for each of the variables x\, . . . , x c i 1 , Z\ , . . . , Zd 2 • Let M be the max- 
imum of the absolute values of these 2d numbers. Thus |x,|, < M for 
all i. Let C be the maximum of the absolute values of all coefficients, and 
r be the number of monomials of f{x). Then 

L := -rCM D < /(x,z) < rCM D =: U , 

as we can bound the absolute value of each monomial of f(x) by CM D . 

(2) Let i := 0. 

(3) Using the algorithm of |Corollary 15[ compute an upper bound u for the 
problem 

max g(x, z) := /(x, z) — L, 
s.t. (x,z) e PD (R dl x Z d2 ) 
that gives a (1 +<5)-approximation to the optimal value. Let t/j+i := L,+m. 
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(4) Likewise, compute an upper bound u for the problem 

max /i(x, z) := Ui — /(x, z) 
s.t. (x,z) e Pn (R dl x Z d2 ) 

that gives a (1 + ^-approximation to the optimal value. Let L i+ i := Ui — u. 

(5) i:=i + l. 

(6) Go to El 

Lemma 17. \Algorithm 16\ is correct. For fixed < S < 1, it computes the bounds 
L n , U n satisfying (|24p and (|25[) m iime polynomial in the input size and n. 



Proof. We have 
and 



Ui-L i+1 < (l + *)(Ui-/(x ininj z mill )) (26) 
< (l + <5)(/(x max ,z max ) -L 4 ). (27) 



This implies 

Ui+i - < m - U) + (1 + <5)(/(x max , Zmax ) - /( 

^min: ^min)) ■ 

Therefore 

U n - L n < 5 n (U - Lo) + (1 + S) S 1 (J(x max , z max ) 

1 - (5™ _1 

= 5"(t7 - Lq) + (1 + 6) _ (/(X roax , Z max ) - /(x m i„, Z min )) 

^ - 7 (/(x max , z max ) /(x m i n , z m j n )) (u > oo). 

The bound on the running time requires a careful analysis. Because in each step 
the result u (a rational number) of the bounding procedure ( jCorollary 15[ ) becomes 
part of the input in the next iteration, the encoding length of the input could grow 
exponentially after only polynomially many steps. However, we will show that the 
encoding length only grows very slowly. 

First we need to remark that the auxiliary objective functions g and h have 
integer coefficients except for the constant term, which may be rational. It turns 
out that the estimates in the proof of lThcorcm 11 (in particular, the local Lipschitz 
constant L and the lower bound on the optimal value) are independent from the 
constant term of the objective function. Therefore, the same approximating grid 
— Z dl x Z d2 can be chosen in all iterations of Algorithm 16 the number m only 



depends on 6, the polytope P, the maximum total degree D, and the coefficients 
of / with the exception of the constant term. 

The construction in the proof of |Corollary 15| obtains the upper bound u by 
multiplying the approximation /(x £ ,z e ) by (1 + 8). Therefore we have 

U l+ i = Li+u 

= L i +(l + <y)(/(x 6 ,z e )-L i ) 

= -5L i + (l + 8)f(x £ ,z £ ). (28) 

Because the solution (x e ,z £ ) lies in the grid ^Z dl x Z d2 , the value /(x c ,z e ) is 
an integer multiple of m~ D . This implies that, because Lq < /(x £ ,z £ ) < Uo, the 
encoding length of the rational number /(x e , z e ) is bounded by a polynomial in the 
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input size of / and P. Therefore the encoding length f/j+i (and likewise Li+i) only 
increases by an additive term that is bounded by a polynomial in the input size of 
/ and P. ' □ 

We are now in the position to prove [Theorem 2l 

Proof of \Theorem 6 A Clearly we can restrict ourselves to polynomials with integer 
coefficients. Let m = (D + l)diA, where A is the number from IThcorcm 81 We 
apply [Algorithm 16| using < 5 < 1 arbitrary to compute bounds U n and L n for 

n= \-log s (2m D (U -L ))]. 

Because n is bounded by a polynomial in the input size and the maximum total 
degree D, this can be done in polynomial time. Now, by the proof of lLcmma 17[ 
we have 

1 - <5" _1 

U n - L n < S n (U - L ) + (1 + S) (/(Xpax, Z max ) - /(x ram ,Z mm )) 

— 2^ ~\ j(/( x max! z niax) /(x m i n > z min)) • (^9) 

If / is constant on P n (R dl x Z d2 ), it is constant on Pn (^Z dl x Z d2 ), then 
U n — L n < \ m ~ D ' ■ Otherwise, bv lTheorem 14i we have U n — L n > /(x max , z max ) — 
/(x m ; n ,z m i n ) > mT D . This settles part (a). 

For part (b), if / is constant on Pn (R dl x Z d2 J, we return an arbitrary solution 
as an optimal solution. Otherwise, we can estimate further: 

U n - L n < Q + JZ^j (/(Xmax, Zmax) - /(x m m, Z min )). (30) 

Now we apply the algorithm of lThcorcm II to the maximization problem of the 
polynomial function /' := / — L n , which is non- negative over the feasible region 
Pn(R dl x Z d2 ). We compute a point (x e <,z e /) where e' = e(± + j^f ) 1 such that 

\f (^c' ) Zg' ) f (x max , Z max ) | ^ € f (x max , Z max ) . 

Then we obtain the estimate 

| fip^-e' i ) /(^maxi 2 max ) | ^ 6 (^(x max , Z max ) -^n) 

<e'(U n -L n ) 

^ ^2 "J fi^J (/( x maxj Z max ) — /(x m ; n , Z m ; n )) 
^(y(Xmax) Z max ) f (x mm , Z max ) ) , 

which proves part (b). □ 
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